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Abstract. We review some modern applications of the theory of few- 
body encounters between binaries and single stars. In particular we focus 
on the treatment of adiabatic encounters, in a regime which is of im¬ 
portance in encounters between a star and a planetary system in a star 
cluster. 


1. Introduction: examples and applications 

Roughly speaking, gravitational scattering will be defined as the study of few- 
body encounters in which particles interact by Newtonian gravity, with certain 
types of initial conditions; namely, a few (normally two or three) bound subsys¬ 
tems, such as single stars or binaries, approach from infinity. Then the problem 
is to characterise the outcome, usually in a statistical sense. 

The new book by Valtonen & Karttunen (2006) will very quickly become the 
standard reference on this problem. While approximate analytical methods yield 
useful results in some limiting situations (see, for example, Sec.2.), computer 
simulation is an essential tool. And while efficient codes exist for computing 
individual scattering encounters (e.g. triple 1 , fewbody 2 ) the scattering packages 
in starlab 3 4 provide the additional, extremely valuable functionality of efficiently 
sampling parameter space. Many possibilities exist for graphic rendering of 
individual encounters, such as GLaninr. 

1.1. Examples and applications 

Though the topic was developed throughout the twentieth century, it remains 
topical, because of fresh applications, such as the following. 

Scattering of normal stars by a binary black hole As a result of mergers, cen¬ 
tral binary black holes are expected in many galaxies. Observationally, binary 
black holes are studied at high energies (e.g. NGC 6240, studied with Chandra: 
Komossa et al 2003) and at visual wavelengths (e.g. the famous Tuorla object 
OJ287: Sillanpaa et al 1988). A problem of long standing is the evolution of 


1 http: //www. ast. cam. ac .uk/~sverre/web/pages/nbody. htm 

2 http: //www. astro .northwestern. edu/~f regeau/code/f ewbody/ 

3 http: //www. ids . ias . edu/~starlab/ 

4 http: //grape . astron. s .u-tokyo .ac.jp/~makino/softwares/GLanim/ 
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the relative orbit of the black holes, as they scatter stars from the surrounding 
galaxy (Hills 1983, Gould 1991b, Mikkola & Valtonen 1992, Fukushige et al 1992, 
Quinlan 1996, Zier & Biermann 2001, Merritt 2001, 2002, Yu & Tremaine 2003, 
Milosavljevic & Merritt 2003, Chatterjee et al 2003, Gualandris et al 2005). 

Evolution of planetary systems in star clusters For the solar system, though 
it is not now in a star cluster, the question of stellar perturbations has been 
considered for a long time (e.g. Lyttleton & Yabushita 1965). In recent years 
this question has arisen because of the initially surprising absence of planets 
(searched for photometrically) in star clusters (e.g. the globular star cluster 47 
Tuc; see Gilliland et al 2000, Weldrake et al 2005). From the theoretical point 
of view the study of such scattering is simplified by the fact that one component 
of the participating binary is effectively massless (Davies & Sigurdsson 2001, 
Spurzem et al 2003, Fregeau et al 2005). We shall say more about this in Sec.2. 

Capture of exotic particles by multiple systems The idea here is that the flux 
caused by interactions with normal matter may be enhanced in situations where 
the particles can be trapped gravitationally, which allows them multiple op¬ 
portunities of interacting. This is a more speculative problem, but has still 
attracted considerable interest (Press & Spergel 1985, Gould 1987, 1988, 1991a, 
Damour Sz Krauss 1998, 1999, Gould & Alarn 2001, Lundberg & Edsjo 2004). 
Theoretically it is the case where the incoming particle is nearly massless. The 
complexity of the problem is illustrated by a comparable problem of solar system 
dynamics: what fraction of comets and asteroids are destroyed by colliding with 
the sun? Several planetary resonances are involved in the fact that almost half 
of a sample of near-Earth asteroids end by colliding with the sun (Farinella et 
al 1994). It is also a common fate of short-period comets (Levison & Duncan 
1994). 

The M4 triple The nearby globular cluster M4 contains a triple system in which 
a distant companion, of mass comparable with or somewhat larger than Jupiter’s 
mass, is in orbit about a binary consisting of a white dwarf and a millisecond 
pulsar. It is likely to have formed in a four-body scattering encounter between 
two binaries (Rasio et al 1995, Ford et al 2000, Fregeau et al 2005). 


Adiabatic encounters 


Consider a three-body scattering event involving a binary, with components 
of mass mi, m?, and a third single star of mass m 3 . In this section we shall 
be concerned with situations where the encounter is both tidal and adiabatic. 
Let the (initial) semi-major axis of the binary be a, and v the typical relative 
velocity of the components. Then we shall estimate v 2 ~ GAR 2 /a, as in a 
circular orbit, where M \2 = m\ + m 2 . Let the relative speed of the star and 
the binary, when far apart, be Vj n f, and let q be the distance of closest approach 
between them. Taking a Keplerian approximation for their relative motion, 

2GM123 

we see that, at close approach, their relative speed is V 2 = V 2 lf -\ -, 

where M 123 = mi + m 2 + m 3 , while the eccentricity of their relative motion is 
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qVj nf _ 
GM\ 23 


By plotting the curves e' 


2 , V/q = v/a (using the foregoing 


estimates) and q = a we distinguish those encounters which are near-parabolic, 
adiabatic and tidal, respectively (Fig.l). 



Figure 1. Regimes of tidal, adiabatic and hyperbolic encounters. En¬ 
counters to the right of the vertical line are tidal. These curves have 
been sketched for the case of equal masses, but do not depend sensi- 
ti vely on th e masses unless m 3 is very large. Vjnf is plotted in units of 
\J GMyxifa. 


Scattering problems like this can be approached analytically in various lim¬ 
iting regimes. High up in the diagram are non-adiabatic, impulsive encounters, 
which covers most situations involving “soft” binaries. At the bottom, just above 
the horizontal axis, in the tidal regime, are near-parabolic, adiabatic encounters. 
A theoretical study of encounters in this regime was carried out by Roy & Had- 
dow (2003), who gave explicit formulae for the change in energy of the binary 
during the encounter. Heggie & Rasio (1996) had done a similar job for the 
change in eccentricity of the binary, and had actually extended their result to 
the regime of hyperbolic, adiabatic encounters. This regime is important for 
one of the problems mentioned earlier (Sec.1.1.), viz. the evolution of planetary 
systems in star clusters, where typical values of Vj n f are comparable with the 
orbital velocity of a planet at a few AU. The main purpose of the present section 
is to extend the principal result of Roy & Haddow to this regime. 

Using first-order perturbation theory, and truncating the perturbing poten¬ 
tial between the binary and m 3 at quadrupole order, Roy & Haddow show that 
the change in energy of the binary is 


be = 


Gm\m2m 3 
Mi-2 


r°° ■ dlZ , 

R.— 


in 


— OO 
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where t is time, R is the position vector of m 3 relative to the barycentre of the 
binary, and 1Z is the perturbing function. In turn, this can be taken to be of the 
form 


R R 5 


|eia 2 (a.R ) 2 - ^e 2 6 2 (b.R ) 2 - i ee 3 a 2 R 2 


cos M + 3 e 4 a 6 a.Rb.RsmM 


(2) 

where a, 6 , e are, respectively, the semi-major axis, the semi-minor axis and 
eccentricity of the binary, a, b are unit vectors along the axes of the orbit of the 
binary, M is the mean anomaly of the binary, and the remaining coefficients are 
defined to be 


ei = J-i(e) - 2eJ 0 (e) + 2eJ 2 (e) - J 3 (e) 
e 2 = J_i(e)-J 3 (e) 
e 3 = eJ_i(e) - 2J 0 (e) + 2J 2 (e) - eJ 3 (e) 
e 4 = J-i(e) - eJo(e) - eJ 2 (e) + J 3 (e), 


in which J n is the Bessel function of the first kind of order n. 

Roy & Haddow proceed to use formulae for parabolic motion for R. We 
follow their same procedure using formulae for hyperbolic motion (e.g. Plummer 
1918), i.e. 


R = a / (e / — cosh F) A + b' sinhFB (3) 

R = a' (e 1 cosh F — 1 ), (4) 

where A, B are unit vectors aligned with the axes of the hyperbolic relative orbit 
of m 3 , a', 6 ', e' are the hyperbolic analogues of a, 6 , e for this orbit, and F is the 
eccentric anomaly, related to time by 


n't = e' sinh F — F, (5) 

in which 

n'V 3 = GM 123 (6) 

and we have assumed t = 0 at the time of closest approach. We also have 
cosAf = ;Rexp(m(t — to)), where to is the time of pericentric passage in the 
binary, and n 2 a 3 = GM\ 2 - 

From eqs.(l) and (2), we see that a typical term in the integrand of 5s is 


R 


d ((a • R) 


3R \ R 5 


cos M = 3? 


2a • Ra • R *R - 5(a • R ) 2 R* 
R6 


F exp (in(t — to)) 


( 7 ) 

where * denotes differentiation with respect to F. Indeed all the terms in eq.(l) 
can be expressed in terms of the integral 


1 = 



f(F)F 

R 6 


exp (in(t — to))dt, 
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Figure 2. The modulus of the integrand of eq.(8) in the region - 
1 < 5iE < 1,0 < ^E < 2 of the complex E-plane, for the case e' — 1. 
Essentially the shape is that of a saddle, interrupted by a pole near 
F = i. 


where / is a polynomial in coshE and sinhE. By means of eq.(5), we then have 


fC 

F = exp(— into) / 

J -! 


m 

R 6 


exp < i 


. n(e' sinh E — E) 


TV 


dF. 


( 8 ) 


We handle this integral in the same way as in Heggie & Rasio (1996). 
Basically the method is steepest descents, but at the location of the saddle of 
the exponent, it turns out that R = 0. (This is responsible for the pole at 
about F = i in Fig.(2).) First, therefore, we clear off the factors of R in the 
denominator. We deform the contour off the real E-axis to avoid the singularity 
at E = 0 in the following equivalent expression: 


T = exo (—into) 


roo 




where we have used the derivative of eq.(4). Integrating by parts we can convert 
R*/R e into 1/(5E 5 ), and in differentiating the remainder of the integrand we can 
treat everything except the exponential as constant: the adiabatic assumption 
implies that n/n' 1, and so the derivative of the exponential factor dominates. 
Thus we have approximately 


E = exp(— into) 


1 in 
5 n'a' 2 e' 


roo f(F) 1 f . n(e' sinh E — E) 

7-oo sinh E R 4 ' ^ \ n’ 


dF, 
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where we have made use of eq.(4) again. Doing this twice more gives the ap¬ 
proximate result 


1 = exp(—into) 


1 

1.3.5 


/ in 
\n'a' 2 e' 



f(F) 
sinh 3 F 


exp 


J .n(e’ sinh F — F) 
n' 


dF. (9) 


Now we deform the contour onto the saddle of the exponent, which occurs 
where its derivative vanishes, i.e. e! cosh F —1 = 0. We choose the root Fq = 
i arccos(l/e / ), which leads to 


sinh Fq = i\jl~ coshN 0 = (10) 

A quadratic approximation of the exponent around F = Fq gives 

i—(e sinh F — F) ~-- ( vV 2 — 1 — arccos(l/e / ) + - Ve /2 — 1(F — No) 2 } . 

n n f 2 J 

Now the integral in eq.(9) is easy, and gives 


1 


exp(—int 0 )——— 
1.3.5 


/ in 
\n'a' 2 e' 



ir i/4 m b) 

sinh 3 Fq 


x exp 



1 — arccos 



x 


Since f(F) is dehned implicitly in terms of R, R and their derivatives (cf. eq.(7)), 
it is helpful to know that, when F = Fq, we have R = a\e' — l/e , )(A + iB), R = 
0, R* = a'l/l — l/e /2 (—iA + B) and R* = ia!\Je' 2 — 1, where we have used 
eqs.(3),(4) and (10), and the fact that b' = a' Ve /2 — 1. (The fact that R = 0 at 
the saddle is the reason for the three integrations by parts.) 

Now we see that the integral (with respect to t ) of eq.(7) is approximately 


K exp (—into) 


iV2n(e' 2 - l) 3 / 4 
15a /3 e /2 



(—5)(a.A + ia.B) 2 


x 


x exp 



1 — arccos 



We use eq.(6) and the corresponding relation for the binary to eliminate n and 
n ', and eliminate a' in favour of the distance of closest approach, q. Thus the 
contribution of this term to 5e in eq.(l) is 


Gmxm^m.zM^ \/2 tt (e! + 1 ) 3 / 4 g 3 / 4 ( c ^ 

10 H 2 ^ 74 (_5ei) 


M- 


5/4 

123 


a. A) 2 — (a.B) 2 ) sinnto— 


—2a.Aa.B cos 



yfe ! 2 


— 1 — arccos(l/e / ) 
(e' - l) 3 / 2 
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For ease of comparison, and a little economy, we introduce K 
so this contribution to 8s is 


2M 12 r/ 3 

Mi 23 a 3 ’ 


and 


Gmim 2 m 3 Jit ie! + l) 3 / 4 ^ /2/ _ 2 » r,,, ; , 2 v^ 2x ■ 

' M 12 g 3 lo 2 3 / 4 e 2 7 (" 5 e i a “) {(( a ' A ) - (a.B) ) smni 0 - 


i. Aa.B cos nto} exp I — 


K Ve' 2 - 1 — arccos(l/e') 
v/2 (e' - l) 3 / 2 


This is identical with the result of the corresponding terms in Roy & Haddow 
(the unnumbered equation below their eq.(17)), except for the factor involving 
e', and the exponent. In the limit e! —> 1 the results are in agreement. 

We can obtain a more convenient final expression by making the same 
changes to eq.(19) in Roy & Haddow, yielding 


Grrnm 2 m :i (e ' + 1 ) 3/4 ^ 5/2 ( I< \Je' 2 - 1 — arccos(l/e') 

AW 8 2 3 / 4 e /2 1 6XP \^2 (e' - l ) 3 / 2 

x jeia 2 [sin( 2 u; + nfo)(cos 2 i — 1) — sin( 2 a; + nto) cos 2 i cos 2 H— 

—3 sin(2a; + nto ) cos 2H — 4 sin 2H cos(2a; + nto) cos i] + 

+e 2 6 2 [sin( 2 a; + nto){l ~ cos 2 i) — sin( 2 a; + nto) cos 2 i cos 2Q — 

—3 sin(2a-’ + nto) cos 2H — 4 cos(2a + nto) sin 2H cos i] + 

+e 4 a 6[—2 cos 2 i cos( 2 a + nto) sin 2 H — 6 cos( 2 a + nto) sin 2 fl — 

—8cos 2 H sin( 2 a; + nto) cos t]} , (11) 



where a, H and i describe the orientation of the path of m 3 in a frame aligned 
with the elliptical orbit of the binary. Thus H is the longitude of the ascending 
node, measured in the plane of the binary from its pericentre in the direction of 
its motion. 

It would be desirable to test this formula numerically. So far the only tests 
to have been carried out involve qualitative comparison of scatter plots produced 
using both this formula and numerical scattering experiments. One reason why 
such a comparison is desirable is the cavalier nature of the perturbation calcu¬ 
lation we have carried out. In particular the end result is exponentially small 
in the ratio a/q , whereas we have neglected terms which are only algebraically 
small (e.g. higher order terms omitted in eq.(2)). Nevertheless the correspond¬ 
ing result for e! = 1 works well (Roy & Haddow 2003). Finally we remark that 
eq.(ll) gives a null result when e = 0. Roy & Haddow show how to deal with 
this when the encounter is parabolic, and it is expected that a result for the 
general case could be provided along the lines of the above calculation. 
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